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Abstract 

We construct a general stochastic process and prove weak convergence results. It 
is scaled in space and through the parameters of its distribution. We show that 
our simplified scaling is equivalent to time scaling used frequently. The process is 
constructed as an integral with respect to a Poisson random measure which governs 
several parameters of trading agents in the context of stock prices. When the trading 
occurs more frequently and in smaller quantities, the limit is a fractional Brownian 
motion. In contrast, a stable Levy motion is obtained if the rate of trading decreases 
while its effect rate increases. 

Key words: fractional Brownian motion, arbitrage, stock price model, stable Levy 
motion, long-range dependence, self-similarity 
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The limit is a fractional 
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Brownian motion (fBm) or a Levy process depending on the particular scaling. While 
the motivation of such analysis originates from data traffic in telecommunications, 
both fBm and Levy processes have recently become prevalent in finance. Thereby, 
we construct a general stochastic process based on a Poisson random measure N, 
interpret it as stock price process and prove weak convergence results. 
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We consider a real valued process of the form 

/OO (•OO poo ' / i \ / \ " 

/ / a n ru /( -/(_ N n (ds,du,dr) (1) 

where / is a deterministic function satisfying Lipschitz condition, (a n ) is a scaling 
sequence, and r and u are marks of the resulting Poisson point process with s denoting 
the time. The mean measure on U is either a probability measure or an infinite 
measure on R + . The process Z n depends on the scaling parameter n £ K+ through 
not only a n , but also through the mean measure fi n of N n which is taken to have 
a regularly varying form in u to comply with the long-range dependence property 
of teletraffic or financial data. After centering the process Z n , we can obtain either 
an ffim or a stable Levy motion depending on the particular scaling of the mean 
measure of N n and the factor a n as n — s> oo. While fBm is a self-similar and long- 
range dependent model, a Levy process has independent increments and self-similarity 
exists without long-range dependence. 

Our main contribution is the generalization of the previous results 
specific linear form for / and 
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31| with a 



29] with an increasing / satisfying some other technical 
conditions, to those with Lipschitz functions. In this case, not only the proofs require 
more work, but we need Lipschitz assumptions on the derivative of / as well. We 
also show that the time scaling used in previous work can be replaced by parameter 
scaling of the distributions of the relevant random variables. The time scaling has 
been interpreted as 'birds-eye' description of a process, which is not necessary when 



the scaling is interpreted in terms of its parameters. Inspired by 



, we unify the 



results for general forms of / with less stringent conditions in some cases. In 
it is noted that a fractional Brownian motion with H > 1/2 can be approximated 
if the pulse / is continuous on K and has a compact support. As an alternative 
extension, we consider continuous / with no compact support while constructing ([1} 
with N n ~ N n — \i n as a centered process. A secondary generalization of previous 
work is the consideration of Z as a more general process than workload which would 
be positive by definition. We allow a signed process through the choice of real valued 
rate r. 

measure is replaced by a general 
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arrival process and a cluster Poisson process, respectively. Ergodicity is required for 
the limit theorems in the general arrival case as well. On the other hand, most of the 
previous studies on scaled input processes are named infinite source Poisson models 
due to the assumption of Poisson arrivals. 

As for application in finance, the process Z n can be interpreted as the price of 
a stock. The interpretation of (JXJ) as a stock price process has been first presented 
in [l|, Our aim is to construct a model involving the behavior of agents that 

can be parameterized and estimated from data, yet having well-known stochastic 
processes as its limits. While the limiting models fit well to financial data, they do 
not involve the physical parameters of the trading agents. Agent based modeling is 
widely used to find a model that best fits stock price processes. In some studies, agents 
are divided into two groups, mostly named as chartists and fundamentalists. In these 
studies, the two agent groups have different demand functions for the stock 14 L 2ll 23 1 
and the price is generally determined via the total excess demand 
In ([1]), the arrival time s, the rate r and the duration u of the effect of an order 
are all governed by the Poisson random measure. Under the assumption of positive 
correlation between the total net demand and the price change, we expect that a buy 
order of an agent increases the price whereas a sell order decreases it. Each order 
has an effect proportional to its volume and duration. The duration of the effect is 
assumed to follow a heavy tailed distribution. This effect starts when the order is 
given, increases to a maximum which is proportional to the total order amount, and 
then starts decreasing until it vanishes after a finite time. Alternatively, its effect 
may last for some time and leave the price at a changed level on and after the time of 
transaction. The logarithm of the stock price is found by aggregating the incremental 
effects of orders placed by all active agents in [0, t}. As a semi-martingale, our process 
does not allow for arbitrage. It is a novel model which is alternative to existing agent 
based constructions that are semi-Markov processes Q . 

We prove weak convergence results for the two different measures used for the 
duration, separately. As for their connection, the infinite measure being the limit 
of the probability measure is proved in Theorem 1. Although a probability measure 
is appropriate for the duration, its limiting form is able to capture the essential 
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tail behavior which represents long-range dependence and self-similarity properties. 
Moreover, the scalings with the limiting measure is simpler to interpret as follows. 
An fBm limit is found as the frequency of trading increases and the price effect of 
the orders decreases as shown in Theorem 3. In contrast, a stable process which is 
a scaled version of a stable Levy motion is obtained in Theorem 5 when the rate 
of trading decreases and the effect of the orders increases. Theorems 2 and 4 are 
concerned with fBm and stable limits, respectively, for the case of probability measure 
for the duration. The hypotheses in these theorems involve extra scaling needed for 
obtaining the limiting measure in addition to the scaling of the rates and effects. The 
stable process obtained in the limit has a skewness parameter that depends on the 
distribution of the rate r as it can take positive or negative values. 

The paper is organized as follows. In Section 2, we define the ingredients of 
the workload model, namely the random variables associated with the process to be 
constructed. The price process is defined in Section 3 with various assumptions on 
the parameters of the Poisson random measure. Section 4 includes limit theorems 
for fractional Brownian motion. Finally, limit theorems for stable Levy motion are 
proved in Section 5. 

2. Infinite Source Poisson Agents 

In this section, we state our main assumptions and notation for constructing a 
stock price process. We assume potentially an infinite pool of agents and gather the 
agents' trading processes under a Poisson random measure. Each agent's trading 
starts according to the underlying Poisson process and it ends after a random amount 
of time. Although the same agents may be returning for further transaction, new 
arrivals are assumed to be independent and identically distributed. 

The agents are called infinite source Poisson as a term borrowed from traffic mod- 
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eling in telecommunications based on Poisson random measures (e.g. [l( 
Since self-similarity and long-range dependence are common statistical properties ob- 
served in both Internet traffic and financial data, the stochastic models also have 
similar features. 
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In [3| , a finite number of identically distributed semi-Markov processes represent- 
ing the trading states of the agents over time are aggregated to form the price process. 
The number of semi-Markov processes, equivalently, the number of agents has been 
taken to infinity only in the limit. On the other hand, our price model is a stationary 
process with an infinite source of agents that arrive according to a Poisson process. 
We concentrate on scaling the other parameters that have physical interpretations for 
obtaining the limiting stochastic processes. 

Let (f2, _F, P) be a probability space. Let Br denote the Borel er-algebra on R. Let 
N be a Poisson random measure on (R x R + x R, Br ® Br + <g> Br) with mean measure 

fi(ds,du,dr) = Xdsv(du)^/{dr) (2) 

where A > is the arrival rate of the underlying Poisson process, 1 < 5 < 2, 7 is 
the distribution of a random variable R and v is either a probability measure that 
satisfies 

u- & 

vidy) ~ h(u) — — as u — > 00 (3) 


where h is a slowly varying function at infinity, that is, h is such that 

lim h(un)/h(n) = 1 (4) 

n— »oo 

or a measure given by 

v(du) = v,- 5 - 1 du . (5) 

Each atom (Sj,Uj, Rj) of N can be interpreted as an order from an agent, where 
Sj is the arrival time of the order, Uj is the duration of its effect on the price, and 
Rj denotes its rate which also plays the role of conversion to monetary units. The 
sign of Rj could be positive or negative depending on the order being a buy or 
sell order, respectively. Under assumption ([5]) for the measure v, the duration U is 
obtained from a diffuse measure on (0, 00) with no probability distribution. This case 
is studied for suppressing the less significant details in the proofs of the convergence 
theorems. In case of ([3]), U follows a heavy tailed distribution with finite mean but 
infinite variance, and the convergence proofs involve the function h. Although the 
latter case is physically more meaningful, the scalings are more involved as well in 
the limit theorems for the self-similar processes fBm and stable Levy motion. 
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Let K : R X M + xl->l denote an effect function such that K(x, •, •) = if 
x < 0. The effect of an order starting at s and ending at s + u depends on the rate 
r of the effect and equals K{t — s, u, r) at time t provided that s < t. The function 
K can also be interpreted as the local dynamics of a transaction as a result of a buy 
or sell order. The rate r will be connected to the quantity of the order which will be 
elaborated further in the sequel. We specify a general effect function / :!->! that 
determines K by 

K{t-s,u,r)=ruf[ t —^-\ (6) 



u 

for t > s, and f(x) = for x < 0. We consider / to be the deterministic shape 
function, or pulse, for the effect which is then shifted to the starting position s, scaled 
and amplified for the duration u, and adjusted once more with the rate/conversion 
factor r, all randomized by N. 

The price process will be constructed as a sum of randomized pulses. Characteri- 
zation of / that will yield an fBm or a Levy motion for the price process is of interest, 
provided that the parameters in ^ are appropriately scaled. In the workload pro- 
, 3, 0, Q> / has the following form 



cesses studied in 
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, x A 1 x > 

m = \ (7) 

x < 

which represents an increasing input, or effect, with unit rate on [0, 1] and remains 
constant thereafter. In [9(, the pulse 

f 1/2- \x - 1/21 < x < 1 

/(*) - I ~ ~ (8) 

I otherwise 

is considered for the aim of approximating an fBm. It has a compact support repre- 
senting a limited effect that vanishes after the duration of the pulse. These special 
pulses, in other words effect functions, are sketched in Fig.l. General Lipschitz func- 
tions are also considered in 

The effect function is left unspecified in [3, H, Q] but with conditions on the 
tail properties of its distribution for large times. The duration is not parameterized 
in contrast to the present work. In [28(, the effect may last indefinitely although it 
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1 



Figure 1: Sample pulses of different kind 



decreases in time with a regularly varying tail. In 22 1. each effect is assumed to 
converge for large times to a finite random variable which has a distribution with a 
regularly varying tail. The stable limits are outlined in [^l with several interesting 
special cases. In 29(, the effect function is deterministic which is randomized through 
a random variable for duration as in our case, but with an additional assumption that 
the effect function itself also has a regularly varyin g ta il. On the other hand, random 
effect functions have been also considered also in 30j where central limit theorems 
are proved under general conditions. As a special case, the effect function could be 
a compound Poisson process as in 
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25|. This could be used to model the buy 



or sell transactions in smaller quantities for a given order in the present work. A 



semimartingale is assumed for the rate of the effect in 



In applications, K can be 



estimated to match the local dynamics of the price change or the workload. 



3. Price Process 



Let the price process be given by Y = {Y(t) : Y(t) = cxpZ(t), t £ R + } where 
Z is the log-price process to be constructed in this section. We aim to introduce 
a stochastic process which is sufficiently general to approximate an fBm or a Levy 
motion, and has an adequate number of physical parameters that can be estimated 
from data. The effect function and the Poisson random measure described in Section 
2 will be the main ingredients. 

Remark 1. Previous models that involve heterogeneous agents usually classify them 



into two separate groups as chartists and fundamentalists according to their trading 
behavior [3J. This can clearly be generalized to several types of agents. The total 
effect Z % from agents of type i, i = 1, . . . ,1, can be further aggregated to form Z as 

n 

z(t) = ^T / z i (t) t>o. 

i=l 

Then, the price Y at time t is given by Y(t) = e z ^ as before. 

We form the log-price process Z by aggregating the randomized effects. More 
precisely, the difference in the effect amplitudes at times and t are integrated with 
respect to the Poisson random measure ./V to yield 



Z(t) = 



oo JO 



(K(t — s, u, r) — K(—s, u, r))N(ds, du, dr) 



t > 0. 



Since the underlying Poisson process has been going on long before time 0, Z has 
stationary increments and Z(Q) — by construction. We think Z(t) — Z(0) as the 
sum of all effects due to all active agents between times and t. We assume that K 
has the form ^ as before, and write 

't - 



Z(t) 



oo JO 



f 



f 



N(ds, du, dr). 



(9) 



The propositions below give the sufficient conditions for Z to be well-defined for a 

finite measure and a specific er-finite measure v in Equation ([2]l. respectively. Let <fi 

denote the characteristic function of Z, that is, </>(£) = Ee^ z for £ G R. 

Proposition 1. Suppose that v is a probability measure satisfying J uv{du) < oo, 
M\R\ < oo and f is a Lipschitz continuous function on R with f(x) = for all x < 
and f(x) — /(l) for all x > 1. Then, Z(t), t > 0, is a finite random variable a.s. 
with characteristic function 



HO = exp 



exp 



i£ru f 



f 



1 > A ds v{du) 7(cZr) 



Proof: The integral J hdfi of a deterministic function h with respect to a Poisson 
random measure defines a finite random variable if J \ h A f | dfi < oo where /i is the 
mean measure. This is clearly satisfied if \ \h\ d/i < oo which also implies that the 



random variable has a finite expectation 
the expression 



Therefore, it is sufficient to show that 



v(du)ds 
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is finite for Z to be well denned as ~E\R\ < oo. Note that 



t - s 



f 



v(du) ds + 



f 



t - s 



v(du) ds 



Considering the two different regions {u : < —s/u < 1 < (t — s)/u} and {u : < 
— s/u < (t — s)/u < 1} for the first integral, and the regions {u : (t — s)/u > 1} and 
{u : (t — s)/u < 1} for the second integral in /, we get 



/(I) - / 



— s 
u 



v(du) 



f 



u |/(1)| v{du) ds 



t - s 



t - s 



v(du) 



-I 
ds 



v(du) 



ds 



Due to the Lipschitz hypothesis on / and since /(0) = 0, we have 



/ < 



l + -)v(du)+ / uM-v(du) 



u M v{du) 



uM- 



t-s 

t - s 
u 



ds 



(10) 



v{du) 



ds 



where M > 0. We apply integration by parts for the inner integrals above. For 
— oo < s < 0, integration by parts yields 



(s + u)v(du) = tF(t - s) 



F{u)du 



and we have t J. v(du) — tF(t — s) where F denotes the cumulative distribution 
function (cdf) of U and F = 1 — F. For < s < t, we get 



Jo 



u v(du) = (t- s)F(t - s) - / F(u)du 



by integration by parts and we have (t — s) v(du) = (t — s)F(t — s). Putting all 
expressions together and using F — 1 — F, we simplify (|T0|) as 



KM 



rt-s 



F{u) duds + M 



J 



F(u) duds 



Changing the order of integration, we get 

/•OO ft — U 

KM 







F(u)dsdu = MtEU 
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which is finite by hypothesis as EU = J uv{du). The characteristic function of Z can 
be found immediately from formulae for integrals with respect to a Poisson random 
measure [26|. □ 

The function / characterized in Proposition [T] represents the local dynamics due 
to the effect of an individual buy or sell order. The change in price, which can be 
non-monotonic, occurs over a finite time U and remains at the same level thereafter. 
The specific shape of / is left unspecified, and so is its sign. In general, we expect a 
buy order to increase the price and a sell order to decrease it. Therefore, if / is chosen 
to be an increasing function, then we could have R > for a buy order and R < for 
a sell order. However, a general form is assumed to leave room for modeling purposes 
in view of real data and to provide mathematical generality. The special case ([7]) used 
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32J,|31| is a linearly increasing pulse as stated in the following corollary. 



Corollary 1. Suppose that v is a probability measure satisfying J uv(du) < oo, 
K\R\ < oo and f(x) = x A 1, x > 0. Then, Z{t), t > 0, is finite a.s. and its 
characteristic function is given by 



exp 



D i£ r(tiA(t-s)+-«A(-s) + ) 



1 



A ds v(du) j(dr) 



for£e 



Remark 2. The log-price process Z which is a semimartingale in general, becomes 
a martingale if its mean is zero. This would be satisfied if Ei? = 0, which corresponds 
to symmetric effects from buy and sell orders, for example. 



The following proposition is based on the results of 
as [0,1] for simplicity without loss of generalization. 



The support of / is chosen 



Proposition 2. Suppose that v(du) — u^ 5 ^ 1 du, Ei? 2 < oo and f : R — > M. is 

Lipschitz continuous with compact support [0,1]. Then, Z is a well-defined random 
variable with characteristic function 



exp 



3 «™[/(V)-/(tt)] - 1 -i^ru [/(*=£)- /(^)]}u -tf-1 Ad»d«7(dr) 



Proof: The Lipschitz assumption on / and that it has compact support [0, 1] imply 
/(0) = /(l) = and 

[f{^)- f{^)fu 1 - 5 duds <oo (11) 
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for each t > 0, and 1 < S < 2 in particular, by |9j, Prop. 3.1]. We sketch the usual 
proof for defining Z as an almost sure limit of zero mean random variables, as in 
Let A k = {2~ k , 2- k+1 ], k = 1, 2, . . . , A = (1, oo) be a partition of R+. Clearly, 

ru f(=?)] u'^Xdsdu^dr) = 0, 

f 

for k = 1,2,..., due to the form of the effect function and boundedness of Ak, and 
also for k — in view of (fTTj) . Hence, the random variables 

r u [/(^)- /(■?)] W(ds, d«, dr) fc = 0, 1, . . . , (12) 

-oo J Ak J — oo 

are well-defined, independent and have zero expectations. One can show that the sum 
of their variances given by 

OO pOO poo 

/ / r 2 u 2 [/(*=•)_ /(^)] 2 u"'- 1 Ads du 7 (dr) (13) 

- oo «/ J — oo 

is finite using (fTTI) and the assumption that Ei? 2 < oo. Therefore, the series 

°° pOO p pOO 

E/ / / ru[f(i^)-f{^)]N{ds,du,dr) 

k—Q J — oo J A k J -oo 

is a.s. convergent by [3, Lemma 3.16]. The limit is denoted by © which is a 
stochastic integral in general since it may be defined by an almost sure limit as above. 
Its characteristic function at £ £ R is the limit of the characteristic functions of 
the partial sums of the random variables (|12p . since almost sure convergence implies 
convergence in distribution. Now, the characteristic function of (|12j) can be written 
as 



exp 



^[/(VHCtt)] _ 1 _ ^ ru [f(t=zy /(_,)] J u-^Xdsdu-yidr) 



since (jl2j) has zero mean. Using the independence of (1121) by disjointness of Ak, 
fc = 0,l,2,..., we can write the characteristic function of their partial sums up to say 
m £ 1i+ as 



exp 



e 'f™[/(V)-/(Tf)] _ l _ l ^ ru [/(*=*)- /(■=*)]} u- s - 1 \dsdu 1 (dr) 



Due to the inequality \e tx — l~ix\ < \x 2 for x £ R and by finiteness of (fT3|) . dominated 
convergence theorem applies and we get the result. □ 
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The pulse © considered in [9'!| is a special case for Proposition[3]where the process 
Z is a zero-mean martingale. In fact, we have E|Z| < oo with the particular pulse 
([5]) which has the shape of an isosceles triangle. In this case, the effect starts at time 
s as the buy or sell order is first given. It increases linearly as the amount traded 
increases, reaches its maximum value at time s + ^ as the highest effect is reached 
and starts decreasing from that point on until it vanishes at time s + u and brings 
the price level back to the original, locally. Such functions / represent the effect of an 
order which takes place over a period, and then vanishes after a while. This scenario 
is a milder and time limited version of the Poisson shot-noise of where a shock 
in the market changes the price through a jump, but then it tends back to its initial 
level by an exponential decay of the first effect. There is no need to get back the 

n 

effect as in the pulses of [9[ for the limit theorems. The effect function may leave the 
price in a level different from the one it found at an arrival, in which case we consider 
the centered process. 

We can show the connection of the two measures in propositions [1] and [5] in the 
next theorem. We introduce a scaling factor n £ Z + which will be taken to infinity 
in the limit. An integral with respect to a probability measure u(du) with regularly 
varying tail converges to an integral with respect to the measure u~ 5 ~ 1 du in the limit. 
We first need the following lemma. 

Lemma 1. Suppose f is a Lipschitz continuous function on K with f(x) = for all 
x < and f(x) = /(l) for all x > 1. Then, 



u 



1 + K 

u~ <5 ~ 1 duds < oo (14) 



'-oo JO 

for 1 < 5 < 3 and n > such that 1 + k > 8. 

Proof: Let us call the integral (p~4|) as /. In view of the assumptions on /, we have 
/ / (u + s) 1+K u- s - 1 duds+ / / u 1+K u- & - 1 duds 

J~ooJ-s JO JO 

t poo rO poo ^ 

/ (t- s) 1+K u- s - 1 duds+ / / t 1+K u- s - l duds\ 

Jt-s J -oo Jt-s J 
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which can be shown along the same lines as in the proof of Proposition [TJ Evaluating 
the above integrals, we find that 

\I\ < M 1+K t 2+K ~ 5 ' ' ' 



(2 + k)(2 + k-5) {2 + k)S (2 + k-8)(1 + k-8) 

1 1 



8(2 + k- 8) 8(8-1) 



□ 



Theorem 1. Suppose that v is a probability measure satisfying J uv(du) < oo and 
has a regularly varying tail as given in the function f : M — > K is Lipschitz 
continuous with f(x) = for all x < 0, f(x) = /(l) for all x > 1 and is also 
differ entiable with f satisfying a Lipschitz condition a.e., and¥J\R\ l+K < oo for some 
< k < 1 with 1 + k > 8 > 1. Let 



[i n (ds, du, dr) 



h(n) 



Xds v n (du) j(dr) 



where v n (du) = v(d(nu)) and 
Z„(t) -- 



oo p OO POO 



oo JO 



u I \ u 



N n (ds, du, dr). 



Then, {Z n (t) — M,Z n (t),t > 0} converges in law to the process 

't - s 



oo r oo poo 



oo JO 



f 



f — 



N'(ds,du,dr) ,t>0 



as n —> oo, where N' = N' — fjf for a Poisson random measure N' with mean measure 
/i'(ds,du,dr) = \u~ s ~ 1 ds du^(dr). 

Proof: For the convergence of finite dimensional distributions of {Z n (t) — MZ n (t), t > 
0}, consider the characteristic function Eexp i^fcLi £,k[Z n (tk) — EZ n (ifc)] for € 
tk > and m 6 N. It is given by 



exp 



oo c oo c oo 



oo JO 



2 ^r =1 ? fc ^[/(^)-/(^)] _j 



fe=l 



/(- 

u 



h(n) 



\dsv n (du)^(dr). (15) 



We first show that the exponent in (fTS"]) is bounded and then use bounded convergence 
theorem to take the limit. This theorem is a generalization of 25J, Thm.l] with the 
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general effect function /. Although we follow the same approach as in 
there are more terms to bound in our case. Let 

't k - s 



25 



Thm.l], 



j(s,u,r) = e 



k ru 



fc=i 



/ 



(16) 

Using the random variable J7, we denote the left hand side of (J3J) as P{[/ > m} below. 
By integration by parts, the exponent in (1151) is equal to 

n s 

d u g(s,u,r)TP{U > nu} \dsdu^{dr) (17) 
h(n) 

where d u is d/du and the hypothesis that v n (du) = v{d(nu)) is used, 
a) Bound for the integrand of (fTT)) for large values of u: 
In view of Potter bounds , for e > there exists % 6N such that 



W{U > nu} 



< 2u~ ,5 max(u" c ,ii e ) 



W{U > n} 

for all n > no and nu > no, that is, u > no/n. Since limn^oo P{?7 > n}n s /h(n) = C 
for some C > 0, we have JP{U > n}n s /h(n) < (C + e) for all n > n' Q for some n' S N. 
Note that C = 1/6 by Assume 

— n o for simplicity of notation. Therefore, we 

get 



P{t7 > ™}4t < 2u" 5 max(w- £ ,w e )(C + e) 



(18) 



for all n > no and u > no/n. 
In ifTfll. we have 



d u g{s,u,r) = i 



'Er=i^«[/(^)-/(^)l 



d u S(s,u,r) 



(19) 



where 



d u S(s 7 u,r) := 



9u 



fe=i 



-/ 



— s 



E^ 



-/ — 

u 



I ( tk — s\ t k 



U \ U I u 



Now, we can bound \d u S\ using the Lipschitz property of / and /' on different regions 
for u and s. Let M > and M' > stand for the Lipschitz constants of / and /', 
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respectively, or their upper bound, whichever is larger. Let us assume M > M 
simplicity of notation. 

i) s < and < s + u < t k 

Since (tk — s)/u > 1 and —s/u < 0, we have 



f 



f 



and 



t ftk — s\ tk- s 



— s \ — s 
u I u 



/(I)-/ 



o-/' 



< M 



1 + 



— s \ — s 
u ) u 



< M 



< M 



due to the form of / and Lipschitz assumptions. Therefore, we get 



> u S(s,u,r)\ < (M 









l+ s - 


+ M 


-1) 


u 




u\J 



;|)E&m = ^Ei&ih 



\d„ 



since |1 + s/u\ < 1 and \s/u\ < 1 in this region, 
ii) s > and s + u < t k 

In this region, /' vanishes and f(—s/u) = 0. Therefore, we have 



\d u S(s,u,r)\ = |/(1)| <Mj2\Zk\ 



hi) < s < tk and t k < s + u 

In this region, f(—s/u) = f'(—s/u) = and we get 
\d u S(s,u,r)\ <2MJ2\&\ h 

k 

iv) s < and tk < s + u 
We have 



tk - 8 



and 



/ ( tk — s\ tk~ s 



U \ U I u 



f — 

u 



—s \ —s 



< M 



< M' 



+ M' 



st k 



The corresponding bound on \d u S(s, u, r)\ follows. 
Now, we can bound the remaining terms in (|19j) by 



2 1 - K El^l K H v |/(^) 



-/ 



—s 
u 



A 2 
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using the inequalities \e ix - 1| < 2 1 " K |x| K and fej \ x j\) < Ej l^'l*' < k < 1 [251, 
and the fact that \e lx — 1| < 2. The index fc is replaced by j in order to distinguish 
the cross products of sums below. We further note that 

tj — s 



f 



/(- 

u 



< (2M) K < 2 K M 



(21) 



since / is bounded and M K < M, assuming M > 1 for simplicity of notation. Putting 
all terms together by (US]), (gDJ), (ED) and i)-iv), we find that (TTT]) is bounded as 

(22) 



\d u g(s, it, r) \ TP{U > nu} A ds du j(dr) 

n /nJ K n ) 



<4M(C + e)^|6 



|r| i?(s, u, tk) max(u e ,u e )u s \dsdu"/(dr) 



where 



B(s,u,t k )= lA^M|^r|rr W K l {s < t , 



L «i, 



tk , |s|*fc 



(23) 



and i?i,fc, . . . , i?4,fc denote the regions in i)-iv). Since [1 A J2j M\£j\ K \r\ K u K l{ s<t .j] 
< 1, we can write 



B(s,u,t k ) < l Rltk +2 



t k . \s\t k 



2- + ^ I l fi4 . fc 



(24) 



lAj2M\^\ K \r\ K u K l {s < t]} l R2 



We keep the extra bounding term for i?2,fe, as the integration in this region is more 
delicate. For fixed k 6 {1, . . . n}, Ri, k , ■ ■ ■ , Ri,k are depicted in Fig. 2. If we choose 
e > such that 

I <8-e <5 <S + e<l + n, (25) 

then the right hand side of (|2"2"j) is finite as as shown in Appendix A. 
b) Bound for the integrand of (fl~7| for small values of u: 

We now consider u < no/n < 1 as n > uq. We use Markov inequality for W{U > 
nu] as in 25], together with the bounds (f2"0)l . ([2~Tj) and i)-iv) in our case. We have 
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Figure 2: Subregions considered for (s,u) 



JP{U > nu} < TEU/(nu) by Markov's inequality. Therefore, we get 

\d u g(s, u, r)\TP{U >n«}^< 2M^^ |r| £ \&\B(s, u, t k ) (26) 

From we can write 

B(s,u,t k )<l Rl . + 2^-^l« 3 . + fe + J^W.* 



j 

considering that u is bounded as u < 1. Now, we have 

1 5+e— 1 1-5-e ^ S+e-1 -£ 1-5 1-S-e ^ 5+e-lj / 1-5 1-5-e / 07 \ 

1 = u u < n Q n n u < n h(n)n u {27) 

since u < uq/ti and n~ e < h(n) for the slowly varying function h when n is sufficiently 



large 
get 



25j . Using ([27)) to increase the righthand side of (l26l and in view of (|24]) , we 



n 4 



|0 uff (s,u,r)|P{£7 > nu} — < 2MEt/n*+^ 1 |r|^" e V |&| 

/i(n) ^ 

V, + 2^-i 1^ + ( 2 £ + 1^*) 1^ + l R2 , k £ Af&l" H" U "l {8 < tj} 

which is integrable over < w < 1, as in part a), in view of the computations in 
Appendix A. 
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As a result of a) and b), the integrand in p7|) is bounded by an integrable function. 
Therefore, we can use dominated convergence theorem to find that 



n s u~ s 
lim TP{U > nu} — - = — 
n-^oo h[n) o 



(28) 



by (J3|) and and then revert (fTT)) by another integration by parts to get the limit 
of ([15]) as 

/OO /*00 /"OO 
/ / g(s,u,r) \u~ s ~ 1 dsdu"f(dr) 
- oo J J — oo 

where g is as in (fTrj|) . It can be shown as in Proposition[3]that the above characteristic 
function and the corresponding process are well defined since \e lx — 1 — ix\ is bounded 
by \x\. Hence, we have shown the convergence of finite dimensional distributions. 

To prove weak convergence in the Skorohod topology on Z?(0, oo), we first observe 
that 



E\Z n (t)-EZ n (t)\ 



l + K 



< 



(29) 



2E|i? 



l + K 



,1 + K 



l + K g 

n 



by 



h(n) 



A ds v n (du) 



25l Lemma 5]. By integration by parts and in view of Potter bounds as before, 



for e > there exists no £ N such that the part of the integral for u > n^/n on the 
right hand side of (|29p is bounded from above by 



/>oo 

2(C + e) / 

J n / 



d_ 

du 



,1 + K 



f 



t - s 



-f 



—s 

u 



\ u - s 



max{u e ,u e } ds du 
(30) 



We have 
d_ 

du 









1 + k" 


u 1+K 













+ (1 + k)u k 



(1 + k)u k 



t- s 
u 

-\ n r 



■f 



n- 

t-s 



II 



l + K 



t - s 



r 



an upper bound for the absolute value of which is given by 



Af 1+K (l + k) {[it _1 (u + s) 1+K + u -1 a(« + s) K ]l fll + u K l R2 
+2u-\t - S ) 1+K l fl3 + u-H 1+K (2 + \s\)l Ri ) 



(31) 
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by Lipschitz assumptions on / and /', where . . . , R4 are as in i) through iv) above, 
with tk = t. Substituting (131 [) in (|30l) and starting the lower limit for u from 0, we 
have an upper bound for the integral in (|30[) given by 

/ [(u + s) 1+K + s(u + s) K ]u- s - 1 m&x(u- e ,u c )duds+ (32) 

-co J-s 

rt rt—a rt k roc 

/ / u K ~ 6 max(u- e ,u e )duds + 2 / / (t - s) 1+K U - s - 1 maxfu" 6 , u e ) du ds 

JO JO Jo Jt-s 

/0 />oo 
/ t 1+K {2 + \s\)u~ S m&x(u~ e ,u € )duds 
-oo J t—s 

which is finite when we choose e as in (|25[) . On the other hand, for < u < no/n, we 
use Markov's inequality as before to get 

n 5 

W{U > nu\ —— < EU n s + e - 1 u- d - c 
h(n) 

Then, the finiteness of the integrals in (|32|) is sufficient again for the integrability of a 
dominating function for < u < no/n < 1 which complements (1301) . It follows from 
dominated convergence theorem that the limit of the right hand side of (f2"9"| exists. 
Therefore, possibly for n > n\ for some n\ £ N, the upper bound in (|29[) is further 
bounded by a multiple of its limit given by 



CiE|i?| 1+K / / u 1+K 







u I \ u 



l + K 



Xu-^duds (33) 



for some C\ > 2. In view of the proof of Lemma [TJ the integral in (1331) is bounded 
by a constant multiple of t 2+K ~ s which clearly dominates E\Z n {t) — EZ n (t)\ 1+K in 
(|29|) for sufficiently large n. Since the increments of {Z n (t) — ~EZ n (t) : t > 0} are 
stationary, this implies that 

E[\Z n {t 2 ) -EZ n (t)\^r \Z n {t) -EZnih)^} < C 2 (t 2 - (t - ti) 2 *^ 

< C 2 (ta - fi) 2+K - 5 (34) 

for < t\ < t < t 2 and some C 2 > 0, by Cauchy-Schwarz inequality and the assump- 
tion that S < 1 + k. This concludes the proof by [a, Thm.13.5 and Eqn.(13.14)] as 
2 + k-S>1. □ 
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Remark 3. Wc consider the scaled measure v n as a scaling of the parameters of 
v. For instance, if v was the Pareto distribution v(du) = Sb 5 ^ 5-1 du for u > 6, 
with parameters 6 > and b > 0, we would have v n (du) = <iu for 

u > 6/n, which would amount to scaling the scale parameter b as b/n. This leads to 
the limiting infinite measure on R + as the cutoff parameter b decreases. Although 
similar interpretations are possible for other probability measures v as well, the scaling 
v n (du) — v{d(nu)) is essentially a time scaling. The random variable U, the duration 
of the sessions, has time interpretation. We look at the workload of individual pulses 
over shorter time periods by scaling U as U/n. 

Note that the limiting process given in Theorem [TJ specifically with / of (J7j); is 
obtained as a limit of sums of scaled renewal processes in [18| as shown in |19j |. 

4. Fractional Brownian Motion Limit 

Fractional Brownian motion is a mean zero Gaussian process Z on R + with Z(0) = 
and covariance 

Cov^r), Z{t 2 )) = ^ (\h\ 2H + \t 2 \ 2H - \h - t 2 \ 2H ) t u t 2 > 



33] 



for ti, <2 > 0, a > and Hurst parameter < H < 1 

In this section, we scale the log-price process as follows to approximate a fractional 
Brownian motion in the limit. Let the rate R be scaled as R/n which can be inter- 
preted as a decrease in the effect of an order in absolute value as n increases. This 
decrease may have arised from an underlying decrease in the volume of the transac- 
tion, for example. On the other hand, we will let the arrival rate A of orders increase 
with a factor which is a function of n. For each n G Z+, let N n denote the Poisson 
random measure with scaled mean measure fi n that involves the scaled arrival rate 
and possibly further scalings. In the following theorems, we prove convergence of 
Z n — EZ„ to fBm with a properly scaled measure v n for a finite measure v as in ([3]) 
and with v as in ([5|) . Compensated Poisson random measure is used when ¥,Z n does 
not exist. 

Lemma 2. Let (x n ) be a strictly positive real sequence and cel. //lim(.x„) = oo. 
then 

c 2 

lim x 2 n (e lc/x " - 1 - ic/x n ) = -— . 

n— J-oo Z 
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— p ic /i 



Proof: Consider the Taylor expansion [17], pg. 184-186] of the function h(x) = e 
around on the disk D = {x G 



c/a 



\x\ < 2}. We have 

■2 



1 - ic/x n ) + 



2x1 



41^2(^)1 



where R2 is the remainder after 3 terms. For fixed r with 1 < r < 2, we can express 
i?2 as 



1 



h(x) 



■ dx 



2111x3 J\ x \ =r x 3 (x — l/x n ) 
By continuity of h, there exists M > such that \h(x)\ < M for all x € D, in 
particular for \x\ — r. Clearly, \x — l/x n \ > \x\ — \l/x n \ = r — l/x n since x n is 
positive. Then, we have the following bound for |i?2(x n )|: 



\Ra(x n )\ < 



Mr 



r - l/x r . 



l/x n 



Taking the limit as n — > 00, we obtain xj; |i? 2 (x„)| 0. 



□ 



Theorem 2. Suppose that v is a probability measure satisfying J uv(du) < 00 and 
has a regularly varying tail as given in the function f : M. — >• M. is Lipschitz 
continuous with f{x) — for all x < 0, f(x) = /(l) for all x > 1 and is also 
differentiate with f satisfying a Lipschitz condition a.e., and ER 2 < 00. Let 



Z n {t) 



00 />00 pOO 



— oc J Q J — 00 



t - s 



f 



—s 
u 



N n (ds, du, dr) 



and 



[i n (ds, du, dr) 



-.2+5 



h(n 



-X ds v n (du) 7(dr) 



where v n (du) — v{d{nu)) and 1 < <5 < 2. Then, the process {Z n (t) — EZ n {t),t > 0} 
converges in law to an fBm with variance parameter 



a z = XER 



OO /-OO 



f 



1 - s 



f 



u 1 s du ds 



as n — > 00. 



Proof: The characteristic function Eexpz^^li £,k[Z n (tk) — EZ n (tk)] for ^ G 
tk > and m 6 N is given by 



exp 



00 roc poo 



1 



-i y\ ik - ■ 

z — J n 



k=l 



tk - S 



/(- 



,2+5 



h(n) 



Xdsv n {du)^(dr). (35) 
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The same approach will be followed as in the proof of Theorem [TJ By integration by 
parts, we find that the exponent of (|35p is given by 

„2+S 



d u g(s, u,r/n) TP{U > nu} , - A ds du j(dr) 

n(n) 



(36) 



Using Potter bounds [(| and Lipschitz conditions on / and /', we get an inequality 
similar to (|22j) for u > no/n given by 

,2+5 

(37) 



d u g(s, u, r/n)\ TP{U > nu} ——— Xds du^(dr) 

n /n J h ( n ) 



<4M 2 (C + e)J2\Ck 



\r\ B(s, u, tk) max(w e ,u e )u s \dsdu"/(dr) 



where B is similar to (|2"3")| but with n = 1 by hypothesis, and e > and % £ N. 
Precisely, 



B(s,u,t k ) = ^1A^M|&| \r\ul {s < tj} 

■ + + 2 ^ U.» + (2£ + lfl 4 , fc 

If we choose e > such that 



K6-e<6<5 + e<2, 

then the right hand side of (|37p is finite along the same lines of the proof of TheorcmQ] 
with k — 1. On the other hand, we can bound (|36p for < u < 1 similarly. Therefore, 
we can use dominated convergence theorem. We have 

„s „,-<5 

lim JP{U > nu} 



as m 



and 



h(n) 



lim n 2 d u g(s,u,r jn) = d u lim n 2 g(s,u,r/n) 

n— >oo n—¥oo 



as g is bounded, hence, uniformly continuous. Then, we get 



lim n 2 g(s, u, r/n) 



(38) 



fc=i 1=1 L v 7 



— s 
it 



/(^) /(- 

u I \ u 
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by Lemma O We now revert (|36p after the limits above, by another integration by 
parts, and get the limit of (f3"5"j) as 



{poo poo poo 2 m m 
-/ / / 2^^^ Cfe7 



/ — 



/(- 



u 2 \dsu d 1 duj(dr) 



which is the characteristic function of Y^TL\ £kZ(tk), where Z = (Z(ti), . . . , Z(t m )) 
is a Gaussian vector with zero mean and covariance 



A :<./,'- 



OO f OO 



u I \ u 



u I \ u 



u 2 ds u 



-S-i 



du. 
(39) 



When (f39l) is evaluated at £ ? - = = 1, the variance coefficient is found to be 



P OO POO 

cr 2 := Var(Z(l)) = AEi? 2 ^ J 



1 - s 



/(- 

u 



ds u 1 s du 



which is finite by Lemma[T]with k = 1. Using the identity 2ab = —(a — b) 2 + a 2 + b 2 
for a, b € R and making several change of variables, we find that the covariance of Z 
in (|39f is given by 



Cov(Zfe), Z(t k )) = °- (if + t 2H - \ tj - t k \ 2H ) 



for tj,tk > with H — (3 — <5)/2. By definition, Z has the characteristic function of 
an fBm. 

Convergence in the Skorohod topology on D(0, oo) follows along the same lines of 
proof of Theorem [T] In this case, we have 



POO rt 

E\Z n (t)-EZ n (t)\ 2 =ER 2 / 

Jo J-, 

and ([M]) holds with « = 1. 



u I \ u 



A ds v n (du) 



As an example, the continuous flow rate model studied in [25] is given by 



Z(t) 



OO p OO />CO 



[{t - s)+ A u - (-s)+ A it] r N(ds, du, dr) 



□ 



(40) 



-oo J J —oo 

with / replaced by the special form (J7J in ([9]). In [25, Thm.2], the limit is studied 
when the speed of time increases in proportion to the intensity of Poisson arrivals. To 
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balance the increasing trading intensity A„ , time is speeded up by a factor n and the 



,5-1 



oo. We can let 



25 



size is normalized by a factor \l/ 2 n( 3 ~ s ' > / 2 provided that X n /i 
A„ = n £+d ~ 1 with e > 0. Taking s = 2, we show the equivalence of the scaling of 
Thm.2] to the scaling in Theorem [5] Note that A„ = n 1+s . The scaled and centered 
process has the form 

Z(nt) - EZ(nt) 



1 

T) 2 



oo poo pOO 



oo JQ 



/ 



nt 



-f 



N n {ds,du,dr) (41) 



where we have written an effect function / in general. Then, we can make change of 
variables s —¥ ns and u — > nu to get 

Zint) - EZ(nt) 



1 



oo JO 
oo poo poo 



oo poo poo 



oo 

r 

— u 

n 



f[ nt-na\ f ,-n» 



-oo JO 

where the mean measure is 



t - s 



nu 



-f[~ 

u 



N n (d(ns) , d{nu) , dr) 
N n (d(ns),d(nu),dr) (42) 



[i n {d{ns), d{nu), dr) 



A„ (nds) h>(d(nu))j(dr) 
n 2+s ds u(d(nu))j(dr) . 



(43) 



In Theorem [2l we start with the scaled process (|42|) essentially. This can be observed 
by the fact that 

N n (d(ns), d{nu), dr) = N' n (ds, du, dr) 
for a Poisson random measure N' with mean measure 



fi' n (ds,du,dr) = [i n (d(ns),d(nu),dr) 



by definition of a Poisson random measure 
scalings in Theorem [2] and 
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Def.V.2.2]. Equivalence of the 



25l Thm.2] is in distributional sense. However, this is 
sufficient for equivalence as the convergence results are in distribution rather than 
almost sure sense. Therefore, we can apply Theorem [2] to obtain the limit as a fBm 



24 



with variance parameter 

poo pOO 

a 2 =ER 2 / [(i _ «)+ a u - Aufdsu- 5 - l du 



(2-S)(3-S) ■ 

It is shown in 25] that the asymptotic behavior of the ratio Xn/n 8 ^ 1 determines 
the type of the limit process when time is speeded up by a factor n. For a choice of 
sequences A„ and n, the random variable §(X n ,n) denotes the number of effects still 
active at time n. It measures the amount of very long pulses that are alive and how 
much they contribute to the total price. The expected value of the random variable 
i(X„,n) is 

Efl(A„,n)~ 



for large n. The limit is considered in the cases where this value tends to a finite 
positive constant, to infinity, or to zero as X n and n go to infinity. We have al- 
ready studied the case of finite constant in Theorem [T] and infinity in Theorem [2J 
the so-called intermediate and fast connection rates, respectively, in view of telecom- 
munication applications. As shown above, our scalings do not involve time scaling. 
They can be physically understood as scalings of the parameters of the log-price pro- 
cess. The slow connection rate will be investigated similarly in terms of the model 
parameters in Theorem [4] in the next section. 

The next theorem is a simpler version of Theorem[2]due to the form of the measure 
v. Note that can be approximated as 

fi n (d(ns), d(nu), dr) ~ n 2+s dsn~ s u~ s ~ 1 du r y{dr) 
= 7t, 2 m _<5_1 ds du"f(dr) 

for large n. This scaling is used below with the simpler form of v. It can be interpreted 
as half way in taking the more involved limit of Theorem [2] 

Theorem 3. Let 



Z n (t) = / / -u 



n 



/(—]-/"* 



u 



N n (ds, du, dr) 



where N = N — ji and 

fj, n (ds,du,dr) = r^Xu^ 8 ^ 1 ds duj(dr). 
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Suppose that Ei? 2 < oo and f : M 
either of the following conditions 



is a Lipschitz continuous function satisfying 



i. f(x) = for all x < and f(x) = /(l) for all x > 1, or 

ii. / has a compact support. 

Then, the process {Z n (t),t > 0}, for 1 < S < 3, converge in law to an fBm with 
variance parameter 



a 2 = AEi? 



OO pOO 



oo JO 



n 1 — )-f(~i 



u 1 s duds 



as n — > oo. 



Proof: Although it can be found from the characteristic function of Z(t) in Propo- 
sition [2] that EZ(t) = for all t > under assumption ii, we form Z as above since 
EZ(i) may not exist with assumption i. For the convergence of finite dimensional dis- 
tributions of {Z n (t),t > 0}, consider the characteristic function Eexp i Y^k=i £kZ n (tk) 
for € R, tk > and m G N. It is given by 



exp 



oo JO 



g(s,u,r /n)n z \u s 1 dsdu"f(dr) 



(44) 



where g is given in (|16p. Note that the characteristic function exists since Z n (tk) are 
well defined in view of (|14[) which follows from Lemma[T]with k = 1 under assumption 
i, and by Proposition [2] under assumption ii. As n — > oo, we will show that the above 
characteristic function converges to 



{/•OO />OC />00 -i m m 
7 / / EEoe 
J-oo JO ./-oo ^ fc=1 =1 



kr 



tk 



u^Xu- 5 - 1 dsdu-f(dr) j(45) 



Due to the inequality |e M — 1 — ix\ < =yx 2 for iGl, the integrand in (|45|) is an upper 



bound to 



\g{s,u,r/n)\n 2 



Therefore, dominated convergence theorem allows us to take the limit inside the 
integral in (|44|) . That is, we must find 



lim e 



£r =1 ^4/(^)-/(^)l 



^ ' T? 



t k - s 



2G 



which is now equal to 

fe=ii=i 



1/ 



7^ (46) 



by Lemma [5J This shows that ([4"4"]l converges to (|4"5|) as n — >• oo by the continuity of 
the exponential function. The variance is evaluated in the proof of Theorem [5] 

To complete the proof, we need to show convergence in D(0, oo) with Skorohod 
topology. This is straight forward since the variance of Z n (t) is already free of n and 
is bounded by a constant multiple of t 3 ~ d by the proof of Lemma [T] □ 

Remark 4. Theorem [3] with condition ii. is Theorem 3.1 of [H where it is noted 
that a fractional Brownian motion with H > 1/2 can be approximated if the pulse 
is continuous and has a compact support. Condition i. above considers an effect 
function which is continuous, but with no compact support as an alternative. 

5. Levy Process Limit 

A process with stationary and independent increments is called a Levy process 



371 ] . The results of this section concerns a particular class of Levy processes, 
namely stable Levy motion [33J. Let Ko = K \ {0}, and let 6 £ (1, 2), and (3 G [—1, 1] 
be the index of stability and skewness parameter, respectively. Then, a 5-stable Levy 
motion L with mean can be defined through its characteristic function 

Ee* L « = exp{~tcr 5 |£| 5 [1 - i/3(sign£) tan(vr^/2)]} 

for (eR, where a > is a scale parameter. 

In this section, we prove that the limiting process is a 5-stable Levy motion under 
different scalings of the price process. Theorem |4] considers a probability measure 
v and Theorem [5] starts with its limiting form. For simplicity of notation, we take 
/(l) = 1 in the effect function. 

Lemma 3. Let N be a Poisson random measure with mean measure 

/.t = Xds u~ s ~ 1 duj(dr) 

and N = N — fi. Then, 



roo roo 



ruN(ds,du,dr) = (A d) 1/s L x {t) + (A C 2 ) 1/S L 2 (t) 
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where C\ — J °° r s ^(dr), C 2 = \r\ s ~/(dr), and L\ and L 2 are independent S-stable 
Levy motions with mean 0, skewness intensity ft equal to 1 and — 1, respectively, and 
scale parameter 

2IY2 - 5) 7n5l 1/5 
■ cos ■ 



5(5-1) 2 



Proof: Putting u' — \r\u, we get 

fOO />oo ft 

r u N(ds, du, dr) 

OO P OO ft 



oo JO JO 



JO JO 
roa rt 



u' [N{ds,d{u'/r),dr) - Xds u'-^r 6 du' j(dr)} 
u' [N{ds, d(u'/\r\),dr) -Xds u'- <5 - 1 |rf du' j(dr)} 
v! [N[{ds,du') - dXdsu'- 5 - 1 du'] (47) 

nt 
v! [N^(ds,du') - C 2 Xdsu'- 5 - 1 du'} 



-m JO JO 

OO ft 



JO 

OO ft 



where we have put 



N{(ds,du') := / N(ds,d(u' /r),dr) 

J{r:r£(0,oo)} 

N' 2 (ds,du') := I N(ds,d(u/\r\),dr) 

J {r:rG(-oo,0)} 



and 

^■oo pO 



poo pO 

C x := / r^idr), C 2 := / 

JO J-o< 



It is easy to verify that both N[ and N 2 are Poisson random measures with means 
CiXds u'^ 5 ^ 1 du' , i = 1,2, respectively. What is more, they are independent as their 
domains are disjoint. Making another change of variable u = u' /{Xd) 1 / 6 for t = 1,2 
in respective integrals in (|4"7|) . we get 

poo />oo pt 

r u N(ds, du, dr) 

/•OO f-t 

(Xd) 1/S / uiNiids^dXd) 1 /^)) -dsu-^du] 
Jo Jo 

nt 
u [N' 2 {ds,d{{XC 2 ) 1/& u)) - dsu- 5 - 1 du] 

--: (XCxf/'L'S) - (XC 2 )^ s L' 2 (t) (48) 



oo JO JO 
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where L[ and L' 2 are independent (5-stable Levy motions with skewness parameter 
P = 1 and scale parameter 

1/5 



for e £ 



B,p 



2r(2 - <5) ttS 

7(^17 C0S T 



g.s 5,156,350]. Now, we have 

i-t 



—L' — Lo :— 



oc JO 



u [N^ds, d(-(\C 2 ) 1/s u)) - ds \u\ 



-s-i 



du] 



^here L 2 is also a (5-stable Levy motion since N^ds.du) :— N' 2 {ds,d{—{\C 2 ) 1 ^ s u)) 



is a Poisson random measure on 
skewness parameter j3 = — 1 
□ 



with mean measure ds I 



\-5-l 



du, but 



33|, pg.5]. We can take L\ = L[ and the result follows. 



Theorem 4. Suppose that v is a probability measure satisfying J uv(du) < 00 and 
has a regularly varying tail as given in the function f : K — >• M is Lipschitz 
continuous with f{x) = for all x < 0, f(x) = /(l) = 1 for all x > 1 and is also 
differ entiable with f satisfying a Lipschitz condition a.e., and E|i?| 1+K < 00 for some 
< k < 1 such that 1 + k > 5 > 1. Lei 



OO />OG />OG 



oo «7 J —00 



n x - al& ru 









-'(?)] 







N n (ds, du, dr) 



and 



[i n {ds, du, dr) 



rrrXds v n (du) j(dr) 



h(n a / s )' 

where v n {du) = v{d{nu)) and < a < S. Then, the process {Z n (t) — KZ n (t),t > 0}, 
for 1 < 5 < 2, converges in law to 

(XER s l {R>0} )^ s L x {t) + (AE| J R| 5 l {i?<0} ) 1 /^ L 2 (t) 

as n —¥ 00, where L\ and L 2 are independent S-stable Levy motions with mean 0, and 
skewness intensity 1 and —1, respectively. 



Proof: The characteristic function of the scaled process is given by 



exp 



OO />00 fOO 



gis^y-^'r)- 



'-00 Jo J -00 Kn a/S ) 
where g is as in (|16l) . By integration by parts, we get 



A ds v n (du) 7(dr) 



00 />oo />oo 



exp 



d u gis,u,n x - a/s r)W{U > nu\ . U ... XdsduMdr) (49) 

h[n a l°) 
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since v n {du) — v(d(nu)). Making a change of variable u to u/n 1 a ' s , the logarithm 
of (gnj) is equal to 

/ / J ^^d u g{ S ,u/n^\n^ s r)W{U >n^ s u} 1 ^ ) Xdsdu 1 (dr) . 

(50) 

Using Potter bounds and Lipschitz conditions on / and /', we get an inequality 
similar to (|22[) . We can bound W{U > n a / s u} n a /h(n a / s ) as in ([T8|) and consider 
\d u g(s, u/n 1 ^ 01 / 6 , n 1 ~ a t s r)\/n 1 ~ a l s separately. As a result, for fixed e > 0, there 
exists no £ N such that for all n with n a / s > uq, we have the following upper bound 
for the absolute value of (|50| when evaluated over u > no/n a / s 



4M(C + e)^2\tk\J J J |r|B'(s,u,i fc ,n)max(u~ e ,M e )u-' 5 AdsdM7(dr) 
where B 1 is analogous to (|23|) satisfying 

B'(s,u,t k ,n)<l Rlikin + J 1 A ^ M\£j\ K \r\ K u K l{ s < t .y j . 1^ (51) 



u / n l-cc/8 «».*.» ^ u/n 1 -"/- 5 tt 2/ n a(i-a/«) ( 

and R\^k,n ■ ■ ■ > R±,k,n are analogous to . . . , i?4,fc with u replaced by u/n 1_Q / 5 . 
The right hand side of (fS"Tj) is integrable when e is chosen as in (|25|) as shown in 
Appendix B. 

For u < no/n a / s , we can find a dominating function for the integrand in (|50|) 
using Markov's inequality. As in the proof of Theorem [TJ we have 

" 2M ^T^T MX) «,**,«) (52) 

where B' satisfies (15T1) . Using (f5Tj) and using an inequality similar to ([27]) in view of 
the assumption u < no/n a ^ s , we can increase the righthand side of (|52[) as 

2AfEf/^+ e - 1 |r|«-'- e £l&l [l fil , fe ,„ + 1k 2 ,,,„ ^M|e,rkr^l{ s < tj } (53) 

fe j 



tk ~ s 
+ 2 u/n l- a /8 1«3 



fn , 1*1** V 1 
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which is integrable over < u < 1 in view of Appendix B. 

We can now use the dominated convergence theorem. Note that 
1 



lim 



d u g{s,u/7i 1 - a / s ,7i 1 - a / s r) = lim d u [g(s, u/n 1 - a ^\ n 1 ^/^)} (54) 



d u lim g{s 1 u/n 



l-a/S „l-a/S. 



where g(s, u/n 1 a / s ,n 1 a / 5 r) is given by 
e ^LL l5 ,™[/(# 7J )-/(^ 7T )]_ 1 



and we have 



lim 

n— f oo 



/ 



f 



' ' -i/n 1 -"/*) " ' 



l{o< s <t fc }(s) 



u j n l- a /S 

(55) 



(56) 



u/n l - a / 5 ) \u/n 1 - a / s / 
To see ([56]). one takes the limit in regions Ri,k, ■ ■ ■ , R±,k, separately. Fig. 3 illustrates 
the function /(•) := /( M / n i~-t/,5 ) over these regions where we consider f(t) — /(0) as 
n — ^ oo. By ([28]l. (|54| . (1551) and <j56j), we take the limit of (J50J) and then revert the 
integration by parts to get the limiting characteristic function as 



exp 



/oo />oo />oo "<■ 
/ / [ e l ^™=i« fcrul {°<=<* fe } - l -i^2^ k rul {0<s<tk} ] Xu-^dsdu-fidr) 
-ocJO J —oo fe— 1 



But, this is the characteristic function of J^fcLi £k%(tk)i where 

/•OO /"OO ^OO 

= III rul {0<s<t }N'(ds,du,dr) 



oo JO 



for a Poisson random measure N' with mean measure /i' = \dsu~ s ~ 1 dwy(dr), and 
N' = N' — //. This characterizes the limiting process by Lemma [3] 

To complete the proof of weak convergence, it is sufficient to show that E,\Z n (t) — 
EZ n (t)\ 1+K < Ct b for some b > 1 and C > in view of the proof of Theorem [TJ In 



the present theorem, we need a finer estimate given in [3a. Lemma 2] and used in 
Lemma 6]. We have 



25 



l\Z n (t)-EZ n {t)\ 



< a 



1-e 



-2/„ 



(57) 



where 
I„ 



1 — cos £ n 



t - S 



fJ, n (ds, du, dr) 
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and a = (J Q (1 — cosx)x 2 K dx) 1 , which is finite with < k < 1. Substituting a n 
and applying integration by parts, we get 

d u k(s, u, n l - a/5 r) V{U > nu} " . A ds 7 (dr) (58) 



h{n a / 5 ) 



where 



/ 



t - s 



k(s, u, r) = 1 — cos ^£ r« 
For latter use, the partial derivative of k in u is found as 

d u k{s, u, r) = sin ( £ur 



/(- 

u 



(59) 



/[ — ]-/ 



u \ u 



f 



f ( t — s\ t — s 



(t) 


— s 






u 





Making a change of variable it to u/n 1 "Z 5 in (|58|) . we get 
1 



7„ = 



■ <9 u fc(s, u/n 1 ^/ 5 , n 1 -"/^) P{C/ > " A ds du j(dr) 



n l-a/S — "V' "/ I ^ ^ J fe^/i) 

Note the similarity of 7 n to (f50|) . Moreover, the inequality sin a; < 2 1 ~ K |x| K A 2 holds 
since since = [e lx — 1 + (e lx + l)]/2 leading to estimates as in (|2ll|) and (l2~Tj) . It follows 
that 

4M(C + e)^JJ J \r\B(s,u,t,n) max^.M 6 )^ Xdsdu-f(dr) 

is an upper bound to |7 n | when it is evaluated over u > no/n a / s where e and no are 
as above and 



B(s,u,t,n) < l Rl ^ +(lAM\C\ K \r\ K u K ).l R2kt . 



1; 



u / n l-a/S 1 u 2/ n 20.-a/S)J R ^ 

For evaluating |7 n | for smaller values of u, we have a bound similar to (|53p . Therefore, 
7 n is bounded by an integrable function uniformly over n in view of the analogous 
computations in Appendix B. By dominated convergence theorem, let 7 = lim„7„. 
We find that 



OO />00 />00 



oo J •* — oo 



[1 — cos(^url{ <s<t})] u 5 1 Xdsdu r y(dr) 
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A 








Figure 3: Graphs of f(x) := f ^ ^l-a/s j m different subregions for (SjU/n 1 a / a ). The function 
/ is arbitrary provided that its derivative and itself are Lipschitz continuous almost everywhere. 



by using the same approach for taking the limit of the characteristic function of the 
finite dimensional distributions above. Then, we can write 

/>oo 

E\Z n (t)-EZ n (t)\ 1+K <a / (l-e- 4/ )r 2 ~ K ^ (60) 

Jo 

for sufficiently large n, by (|57[) . since f — e~ x is increasing in x. Simplifying / further, 
we have 

/OO POO POO 
/ [f -cos^ur^u-^du-yidr) = \t(, 5 E\R\ 5 / (f - cos u) vT^du 
-oo JO JO 

where the second equality follows by a change of variable u to u/(£|r|). Define the 
constant C so that I =: Ct^ s . Now, substituting / in (|60|) and changing £ to , 
we get 

POO 

E\Z n (t) -EZ n {t)\ 1+K < at 1 ^ J (l - e- A ^) C 2 ~ K d( 
which concludes the proof as (1 + k)/5 > 1. □ 

Remark 5. Note that the stable process obtained in the limit is stable with a skew- 
ness parameter that depends on the distribution of the rate R. Moreover, it has 
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stationary and independent increments. Therefore, it is also a ^-stable Levy motion 



33l . Def.7.5.1], but with scale parameter 

aX 1 / 5 (C 1 + C 2 ) 1/5 

and skewness parameter 



Ci + C 2 

by [H pg.s 10,11], where d = ER S 1 {R>0} and C 2 = E^l^o}, see also @, pg.217]. 
In the context of supply and demand, one can interpret C\ as the skewness caused by 
demand and C 2 by supply since they are expected to increase and decrease the price, 
respectively. 

Remark 6. The weak convergence result given in Theorem [4] is proved with Skoro- 
hod's J\ topology. In 25], the analogous result based o n (TQ has been omitted. The 
convergence is shown with M\ topology instead of J\ in [361 ] where the effect function 
is assumed to be monotone increasing in the context of workload input to the sys- 
tem. The authors heuristically argue that some of the individual loads are too large 



36l . Rmk.4.2]. On the other hand, ^22] proves weak convergence with M\ topology 
considering that the limit process has jumps. However, J\ topology also works as 
shown above. The interplay between J\ and M\ is discussed in [2j for sums of moving 
averages. It is proved that J\ convergence cannot hold because adjacent jumps of this 
process can coalesce in the limit. An intuitive explanation is given as the jump of the 
limiting process occurring from a staircase of several jumps. Under certain conditions, 
Mi convergence is shown instead. We have a simpler situation where each arrival of 
the scaled process generates a jump of the limit process as evident from (|56|) . 



The following theorem is based on the simpler form of the mean measure. 



Theorem 5. Suppose the function f : M — > M is Lipschitz continuous with f(x) = 
for all x < 0, f(x) = /(l) = 1 for all x > 1 and is also differ entiable with f satisfying 
a Lipschitz condition a.e., and E|i?| 1+K < 00 for some < K < 1 with 1 + k > 6. Let 



Z n (t) = 



OO /"CO />CO 



/ 



t - s 



-/(- 

u 



N n (ds, du, dr) 



-00 Jo 

where N = N — ji and 

/i„(ds, du, dr) — A n~ s u~ s ~ 1 ds du"/(dr). 

Then, the process {Z n (t),t > 0}, for 1 < 6 < 2, converges in law to 

(XER 5 l {R>0} )^ 5 L 1 (t) + (AElirf l^oj) 1 / 5 L 2 (t) 

as n — > 00, where L\ and L 2 are independent S-stable Levy motions with mean 0, and 
skewness intensity 1 and —1, respectively. 



Proof: Wc will give only a sketch of the proof due to its similarities with the previous 
theorem. The characteristic function for the finite dimensional distributions of Z can 
be written as 

/OO />OC />oo 
/ / g{s,u,nr) An~ d u~ s ~ 1 ds du r y(dr) 
-oo JO J —oo 

with g as in (fT6]l . Making a change of variable u to u/n, we get 

/OO /"CO /* oo 
/ / g(s,u/n,nr) \u~ 5 ~ 1 dsdu"/(dr) (61) 
-oo «/ J —oo 

Now, g(s,u/n,nr) is similar to (|55[) and we take a similar limit to (|56[) with n 1 1 ' 6 
replaced by n. This is justified by dominated convergence theorem since the integrand 
in (|6ip can be bounded as in the proof of Theorem^ Convergence in D(0, oo) follows 
along the same lines, this time with k(s,u/n,nr) in k of (I59I) . □ 

The simpler form of scalings in Theorems 3 and 5 facilitate neat interpretations 
in terms of the parameters of the price process. In Theorem 3, A is scaled as n 2 X 
and r is scaled as r/n, which means that the trading occurs more frequently, but in 
smaller quantities and yields a fractional Brownian motion limit. In contrast, a stable 
process is obtained if the rate of trading decreases while its effect rate increases since 
A is scaled as X/n s and r is scaled as nr in Theorem 5. 
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Appendix A 

We show that the right hand side of (|2"21 is finite. When the right hand side 
of (f2"2")l is splitted over different regions, checking the finiteness of the integrals over 
Ri,k,R^,k,Ri,k reduces to showing that 

rtk—s rt k roo 

u~ s max(u~ e ,u e )duds + / / (tk — s)^ 6 " 1 max(u~ c , u c ) du ds 

Jo Jt k -s 
M poo 



l [ (~ + ^4V" <max ( u e ,u e )duds 

J-oo Jt k -s \ u u J 



is finite. This is indeed true when we choose e > such that 

l<S-e<8 <8 + e<2. (62) 

In region i?2,/c, we have 

/:=// / |r|[l A^2M\^\ K \r\ K u K l {s < tj} ]u- s mnx(u- e 7 u e )dudsj(dr) (63) 
J Jo Jo . 3 

If tj > tk, it can be observed from Fig. 2 that the integral reduces to that over region 
R2.k- If tj < tk, then the integral over i?2,fc yields an upper bound. That is, we can 
replace l{ s <t } by the constant function 1 and get 

I < 1E\R\ 1+K J2 I6T / / L u K ~ s m&x(u- e ,u e )dsdu (64) 
j Jo Jo 

pt k rtk 

+ TE\R\ / u~ s ma,x(u- f ,u e )dsdu 

Ju Jo 

where u denotes a cutoff value of u such that J2j ICjl^^ 1S too large in (|64|) . and we 
use the fact that tk — u < tk for u > after changing the order of integration for u 
and s in (|53"]) . Then, the right hand side of (fM]) is finite if we choose e > such that 

l<6-e<5<5+e<l+K 

which clearly satisfies (1521 since k < 1. 
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Appendix B 

In this part, we show that (|51j) is integrable with respect to max(«" e , u e )u~ s du ds. 
Substituting the limits of integration in regions Ri.k.n, Rs,k,m Ri,k.n shown by I\, I3, 1^, 
respectively, we have 

/0 r-htk—ns 
/ max(«~ e , u e ) u~ s du ds 

-co J — hs 

(htk — ns) max(M _e , u e ) u~ s ~ x du ds 

,tk —hs 
/-oo 



L 



4 



f° f°° ( 2 Isl \ 

/ / ( — ; — I — I Tn&x(u~ e ,u e )u~ s duds 

J-ooJnt k -hs\u/n u^/n^J 



where we put n = n 1 a / 5 . The integrals Ji, J3, I4 are finite for 

l<<5-6<<5<5 + 6<2 



since 

/ u - s duds = C 17 |— < C x t\7 & 

(nt fe - ns)u- & - x duds = C 2 -^— < C 2 t^ 5 

.A k -hs fl 5 - 1 

-ocJhtu-nsKu/n u 2 /n^J n<5-i 

< c3(*r 5 +*^+4^) 

for 1 < (5 < 2 and n > 1, where Ci, C2, C3 £ K. In i?2,/c,n, we have 

nnt-iis I \ 
I 1 A M\^\ K \r\ K u K l {s < t]} J max(u~ e ,u e )u~ s duds 

As in Appendix A, we consider two intervals [0, u] and (u, tk] to evaluate this integral. 
Over the first interval, it is finite for 1 < 6 < 1 + k, and over the latter, it is 
proportional to n 1 ^ 5 which is bounded by 1. As a result, (|5Tj) is finite if we choose 
e > such that 

l<5-e<5<5+e<l+n. 
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